On the asymptotic expansion of the correlators in the XX spin chain. 
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Abstract 

Using the XX- spin chain as an example we show that in general the asymptotic 
expansion of the correlator at large distances is not given by the sum over the low- 
lying particle-hole intermediate states. Only the first two terms in the expansion 
agree with the predictions of the Luttinger liquid theory. The other terms are 
in general given by the intermediate states with the particles and holes with the 
quantum numbers far away from the Fermi- points (at the distances of order of 
| the length of the chain). We argue that the whole expansion cannot be described 

by the Luttinger liquid theory. 



1 1. Introduction 

>' 

in 

Calculation of correlation functions in the ID quantum liquids or spin systems re- 
mains important problem both from theoretical and experimental points of view. Al- 
though the values of the critical exponents corresponding to the power-law decay at large 
distances obtained with the help of the mapping to the Luttinger model (bosonization) 
[1],[2] or conformal field theory [3], [4] are known for a long time, the calculation of the 
constants before the asymptotics (prefactors) remains an open problem. Recently some 
^ ■ progress in this field was achieved in relating of the prefactors to the certain (lowest) 
formfactors of local operators by means of of the conformal field theory [5] . Using the 
Luttinger liquid theory (bosonization) the same results where obtained in [6] . Moreover 
the universal form of the particle-hole formfactors of local operators for the low-lying 
states was obtained [5]. The results [5] are confirmed by an explicit calculation of the 
leading order term of the correlator for the XXZ- spin chain in the magnetic field [7] 
and an explicit calculations of the formfactors for the low-lying particle-hole excitations 
[8] for the XXZ- spin chain in the magnetic field. 

In the present paper we study different formfactors for the XX- spin chain. Our 
results are based on the expression for the general formfactor for the XX- spin chain 
obtained in Ref.[9] in the form of the product depending on the momenta of the eigen- 
states similar to the Cauchy determinant. We calculate both the lowest formfactors 
corresponding to the terms of an arbitrary order in the asymptotics of the correlator 
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and the particle-hole formfactors for the low-lying states corresponding to an arbitrary 
lowest formfactor. In each order we find the particle-hole formfactors in agreement with 
the predictions of the Luttinger liquid theory [5], [6]. Comparing the asymptotics of the 
correlator predicted by the Luttinger liquid theory with the exact results we argue that 
the intermediate states with the particles and holes far from the Fermi- points also give 
the power-law terms in the asymptotics of the correlator. Only the first two terms in 
the asymptotics are given by the low-energy excited states and can be described by the 
Luttinger liquid theory. 

In Section 2 we briefly review the theory of the Luttinger liquid and present the 
derivation of the results [5] on the universal relations for different formfactors in the 
XXZ- spin chain. In Section 3 we calculate the lowest formfactors corresponding to the 
terms of an arbitrary order in the asymptotics of the correlator in the XX- spin chain. 
The expression for the correlator obtained in this section is equivalent to the contribution 
of the low-lying intermediate excited states described by the Luttinger liquid theory. In 
Section 4 we calculate the particle-hole formfactors for different terms in agreement 
with the predictions of the Luttinger liquid theory. Finally in Section 5 we compare 
the predictions of the Luttinger liquid theory for the correlator with the exact results 
obtained by the rigorous method. 

2. Universal relations for the prefactors of the correlator. 

To explain the method of the calculations let us briefly review the derivation of 
the scaling relations for the lowest formfactors and the expressions for the particle-hole 
formfactors in various ID models [5], [6]. For definiteness as an example we will take 
the general XXZ- spin chain in the critical region with the Hamiltonian: 

H = J2 0,X + i + °?<&i + A^K+i) , (1) 
i=i 

where the sites L + l and 1 are coincide. It is well known (for example see [2]) that the 
low-energy effective theory for this Hamiltonian is given by the Luttinger liquid theory. 
In this model one usually defines the fields Ni j2 (x) corresponding to the excitations 
around the right and the left Fermi-points: 

JVi, 2 (x) = j E ^Me-^e-^l/ 2 , Pl (p) = £ ai +p a k , p 2 ( P ) = £ 4 +p c k , (2) 
*-* P ^o P k k 

where ak(ck) are the quasiparticle fermionic operators corresponding to the right (left) 
branch of the Luttinger model and a — > is the regularization parameter. The critical 
exponents of the correlators and other universal relations are determined by a single 
parameter £ which is defined through the ground state energy in the sector with the 
total number of particles and the momentum AN = ANi + AiV 2 , AQ = AN\ — AiV 2 , 
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where AJVi^ - are the numbers of additional particles at the two Fermi points, according 
to the relation 



TT 



AE = -v 



((AN) 2 + (l/fl(AQ) : 



where v -is the speed of sound. For the XXZ- spin chain (1) the parameter £ equals 
£ = 2(7T — i])/tt, where A = cos(?7) [10]. The general operator equation for the operator 
a~ of the original model have the following form: 

a - = ^c' K 1 (K+K 2 ) m e~ i2pFmx e iw ^^ 1 ^~^ 2(x ^ (3) 

m 

where C' m are some non-universal constants and Kf 2 are the Klein factors - the operators 
which commute with the operators Pipip) and create the single particle at the right (left) 
Fermi- point if they act on the ground state. The general operator equation (3) have 
the following sense. The matrix elements of the operator a~ for the low-lying states 
in the original model (1) coincide with the matrix elements of the right-hand side of 
the equation (3) for the corresponding eigenstates in the Luttinger model. If at large 
distances the correlator G(x) = (ct^o-q) is determined by the low-energy intermediate 
states, it is exactly equal to the correlator of the operators at the right-hand side of 
eq.(3) in the Luttinger liquid theory. In particular, according to this prescription, let us 
consider the formfactor of the operator a~ between the ground state \t) (the eigenstate 
with M = L/2 particles (up-spins)) and the eigenstate |A(m)) (the eigenstate which is 
obtained from the ground state with M — 1 particles | A) by moving m particles from the 
left to the right Fermi- point). The matrix element at the right-hand side of eq.(3) for 
the Luttinger model can be easily calculated and we obtain: 

C' m = (\(m)\a;\t)(L/2naf^/Z. (4) 

Calculating the correlator G(x) with the help of the equation (3) and taking into account 
the equation (4) we obtain the general expression for the correlator 

r>( \ l + -\ r< cos(2p F mx) 

CM = W > = E C (M „ /L))(/ ^ M . <«> 

where the Fermi- momentum p F = n/2 and the prefactors C m satisfy the following scaling 
relations for the lowest formfactors (see [5], [6]): 

(-1)-C m / 2 \5/2+- 2 (2/C) 



where (A(m)| is the eigenstate obtained from the ground state (A| by creating m extra 
particles at the right Fermi- point and removing m particles from the left Fermi-point. 
The scaling relations (6) are used to obtain the prefactors of the correlator for the XX- 
spin chain in Section 3. We refer to the equations (3), (5) as to the predictions of the 
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Luttinger liquid theory. In fact as we shall see later the equations (3)-(6) are equivalent 
to taking into account the low-energy excited intermediate states in the correlator. One 
can argue that the contribution of the irrelevant operators to the equation (3) which 
also give the power-law terms in the correlator is equivalent to the contribution of the 
intermediate states with the particles and holes far away from the Fermi-points (high 
energy excited states). We show that the expression (5) does not agree with the exact 
result. Only the first two terms in this expansion agree with the the first two terms in 
the exact expansion of the correlator. That means that in general the low-energy excited 
states do not give the whole asymptotic expansion of the correlator. The main result 
of the present paper is that only the first two terms in the expansion correspond to the 
low-lying particle-hole excitations i.e. can be described by the Luttinger liquid theory. 

Now let us present the derivation of the expressions for the particle-hole formfactors 
i.e. for the formfactors corresponding to the particle-hole intermediate states [5], [6]. 
Below we present the calculations for the particle-hole formfactors corresponding to 
the leading term in the asymptotics of the correlator. The particle-hole formfactors 
corresponding to the other terms in eq.(5) can be obtained in the same way starting 
from the equation (3). Let us denote by (X(pi,qi)\ the eigenstate obtained from the 
ground state (A| by creating the holes with the momenta qi and the particles with the 
momenta Pi with respect to the Fermi-momenta (i — 1, . . .n) located in the vicinity of 
the right Fermi- point. Calculating the matrix element at the right-hand side of eq.(3) 
we obtain: 

„ 2ir Pl(p) /— 

(HPi,Qi)Wo\t) = C{ Pl , qi \e a -U>o — 10>, a = (7) 

where C is the value of the lowest formfactor C = (A | cr^ 1 1) and {pi, qi\ is the eigenstate of 
the Luttinger model with n particles with the momenta Pi with respect to the Fermi-point 
and n holes with the momenta with respect to the Fermi-point at the right branch 
of the Luttinger model. The matrix element at the right- hand side of the equation (7) 
was calculated in Ref.fll] (see Appendix): 

<p,,g,|e a ^>o^|0> = F a (p t ,q t ) = det ij (— ^— ) f[ f + (Pi) f[ /"(ft), (8) 

\Pi Qj J i=\ i=i 

where 

np) = 3£±4, r(9) - ni-i-a) 



r( P )r(o.)' r(i- 9 )r(i-o)' 

In the equation (8) Pi and % are assumed to be integers (corresponding to the momenta 
2irpi/L and 27rgj/L) and Pi > 0, qi < 0. Thus we have calculated the particle-hole 
formfactor in the form (X(pi,qi)\cro\t) = CF a (pi,qi), and have shown that the constant 
C in front of the function F a (pi,qi) is the lowest formfactor. To sum up the formfactor 
expansion for the correlators in the framework of the approach of Ref.[8] (where the 
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dependence on the particle and hole positions was found for the formfactors of the XXZ- 
spin chain in the magnetic field) it is important to have the formula for the sum: 



E E l^,ft)lv^)WL = (9) 



n Pi>0,q;<0 



where p = I]" =1 Pi, <? = Yh=i Qi anc ^ P») 9» are integers. In the same way one can consider 
the particle- hole excitations in the vicinity of the left Fermi- point. Calculating the 
corresponding matrix element at the right-hand side of eq.(3) we obtain: 

„2jt P2(p) l~ 

(HPi,Qi)Wo\t) = C( Pi ,q t \e c -^ P <° — \0), c=^/2, (10) 

where C is again the value of the lowest formfactor C = (A|cr^"|t). The matrix element 
at the right-hand side of Eq.(10) is calculated in the Appendix: 

^— n/ + (Pi)n/"(9<), (ii) 

Pi QjJ i=l i=l 

where now 

1 [P) r(- P )r(i-c)' 1 {q) r(i + g )r(c)" 

In the equation (11) pi and qi are assumed to be integers (corresponding to the momenta 
2npi/L and 2nqi/L) and pi < 0, qi > 0. To sum up the formfactor expansion for the 
correlator it is important to have the formula for the sum: 

E E \FcM\v-^" L = a _ e _* (i2) 

where p = Y^l=\Pii Q — SILi % an d Pi> 9* are integers. The total formfactor is a product 
of the two factors F a {p i) qi) and F c (pi, qi) corresponding to the right and the left Fermi- 
points. The particle- hole formfactors corresponding to the higher order terms in eq.(5) 
can be obtained in the same way starting from the equation (3). 

3. Formfactor for the XX- spin chain and the calculation of the correlators. 

Let us consider the particular case of the XX- spin chain (A = 0): 

H = E 0*X + i + ^f + i) , (13) 

2 = 1 

where the sites L + l and 1 are coincide. The solution of the model (13) has the following 
form [12], [9]. We assume for simplicity L- to be even and M = L/2 to be odd (S z = 
for the ground state and (M — 1)- even, we also assume L to be even so that the ground 
state is not degenerate). Then each eigenstate in the sector with M particles (up-spins) 
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is characterized by the set of the momenta {p} = {pi, . . -Pm} such that pi = 2irni/L, 
rii G Z and each eigenstate in the sector with M — 1 particles is characterized by the set 
of the momenta {q} = {q±, . . . qM-i}, q% = 2n(ni + 1/2)/ L, rii G Z. The ground-state in 
the sector with M particles (up-spins) is given by the configuration {p} = {p ± , . . -Pm}, 
Pi = 2n/L(i — (M + l)/2), (M- is odd), and the ground state in the sector with M — 1 
particles {q } = {qi, ■ ■ ■ qM-i}, qf^ = 2-njLii — M/2). Equivalently one can take the 
shifted momenta 

Pi = (2n/L)(i), i = l,...M, qf> = (2tt/L)(j + 1/2), j = l,...M-l. 

Thus the ground states in the sector S z = (M = L/2) corresponds to the momenta 
(27r/L)n while the eigenstates in the sector S z = —1 (M = L/2 — 1) corresponds to the 
momenta (2ir / L)(n + 1/2), where n is an integer. In terms of the sets of the momenta 
{p} and {q} the formfactor ip({q\) = {{^WoUp}) can be represented in the following 
form [9]: 



nWJ v/IVJ V 6 J n M sin(fe-g,)/2) • U4j 

Let us note that the expression (14) for the formfactor is valid both for the ground state 
configuration of the momenta {q} and for an arbitrary excited state characterized by 
the momenta q±, ■ ■ ■ qu-i- The advantage of the expression (14) for the formfactors in 
comparison with the determinant expression is that with the help of this expression one 
can calculate the formfactor as a function of the set of the momenta {q} for the low- 
energy particle-hole excited states and then calculate the asymptotics of the correlator. 

To calculate the formfactor corresponding to the m-th term in the expansion of 
the correlator (5) we should calculate the expression (14) for the set of the momenta 
{p} corresponding to the ground state and the set of the momenta {q} shifted to the 
right by m in the units 2n/L with respect to the set {q} corresponding to the ground 
state. Clearly, only the denominator gets modified in the process of this shift. Using 
the fact that only the particles close to the Fermi-points change their quantum numbers 
(m « L) and for small arguments the sinuses can be replaced by their arguments the 
resulting formfactor can be easily calculated: 

m m 

Tpm = ik II ((l/2)(3/2) . . . (k - 1 + l/2)(n/L) k ) TJ ((l/2)(3/2) ...(*- 2 + l/2)( 7 r/L) fc - 1 ) , 

k=l k=2 

where ipo is the value of the formfactor corresponding to the ground state (in the nota- 
tions of Section 3 ip = (A|oo"|i), ip m = (A(ra)|<7^"|i)). Then, calculating the products we 
obtain: 

(F(k + l/2)\ 2 T(m+ 1/2) 



D "V r(i/2) J r(i/2) 
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where ipo is again the value of the lowest formfactor. The product in this equation can be 
expressed through the Barnes G-function defined by the relations G(l + z) — F(z)G(z), 
G(l) = 1. The formfactor takes the form: 

^ = Hl) [ 0FG(l/2) ) (,r)"-V2 • (15) 
Using the expression (15) the prefactors C m (5) for m > can be calculated: 

According to the formulas of Section 3 the square of the formfactor |t/> | 2 can be expressed 
through the prefactor of the leading asymptotics as |^o| 2 — (2/ 'L) 1 / 2 Cq, C /2y/n = 
0.147088... [9]. Taking into account the explicit expression for the leading prefactor Co, 
C = 2~ 1 / 2 7r 3 / 2 (G(l/2)) 4 we find that the general prefactor takes the form: 

c m = (-ir^ir^ 2m2 - 2m+1/2 (°( m + l ' 2 ^ ^ m + l i 2 ^ 2 > m > °- 

The equation (17) is the central result of the present paper. In thermodynamic limit 
the asymptotic expansion of the correlator takes the following form: 

C ( 1 ^ cos(7rma;)\ 

G{x) = 7^ + £ ym ^7^ J ' (18) 

where the coefficients y m are equal 

V ( ir 1 ( G(m + l/2) V iY(m + l/2)f 

The first coefficients are: y 1 = —1/8, y<i = 9/2 15 . For small values of m y m decreases 
very fast, however for large m they begin to increase very fast due to the behaviour 
of the function G(m + 1/2). We did not have the explanation of this growth. As we 
will show in Section 4 the results (17)- (19) are equivalent to the contribution of the 
low-lying intermediate excited states to the correlator which will be compared with the 
exact results for the first few terms in Section 5. Let us note that the same expression 
could be obtained with the help of the generalized Fisher-Hartwig conjecture [13], [14]. 

4. Calculation of the particle-hole formfactors for the XX- spin chain. 

Let us calculate the particle- hole formfactors for the state |A(m)) (the set of the 
momenta {q}) corresponding to m extra particles at the right Fermi-point and m extra 
holes at the left Fermi-point. Let us start with the calculation of the formfactor with 
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one particle and one hole at the right Fermi-point. Let us denote the positions of the 
particle and the hole with respect to the first filled level by the positive integers p and 
h (p > 0, h > 0) in the units 2tc/L. The creation of the particle- hole pair leads to the 
modification of both the numerator and the denominator of the equation (14). 

Substituting the corresponding set of the momenta {q} into the expression (14), 
taking into account the condition p, h « L, using the property sin((7r/2 + n/2)/2) = 
1, we obtain after the cancellation of the identical terms in the numerator and the 
denominator the expression 

_ (l/2)(3/2) . . . (p + m - 1 - l/2)(l/2)(3/2) ...(h-m + 1/2) 

iPm{P,h)-y m (p + h)h\(p-l)\ 

where ip m is the value of the lowest formfactor considered in the previous section. Note 
that all factors (n/L) are cancelled. We considered the case h > m. One can prove that 
the case h <m leads to the same expression. Introducing the hole momenta q = — h < 0, 
the last formula takes the following form: 

1 T(p + a) r(l — q — a) . . 

(p - q) T(p)T(a) T(l - q)T(l - a) 

where the parameter a equals 

1 

a = h m. 

In the same way one can easily prove that for n particles and n holes the formfactor 
takes the form ip m F a (pi,qi) where the function F a (pi,qi) is given by the equation (8) 
with the same parameter a = —1/2 + m. Thus we observe that the expression for the 
total formfactor coincides with the prediction of the Luttinger liquid theory presented 
in Section 3 according to the operator equation (3). In fact the equation (3) gives 
a = — v^/2 + m/^ which coincides with our value of a at £ = 1. At m = the 
formfactor corresponds to the leading term in the asymptotics of the correlator. 

Let us proceed with the calculation of the formfactor with one particle and one hole 
at the left Fermi-point. Again we denote the positions of the particle and the hole with 
respect to the first filled level by p and h (p > 0, h > 0) in the units 2n/L. Substituting 
the corresponding set of the momenta {q} into the expression (14), taking into account 
the condition p, h « L, we obtain after the cancellation of the identical terms the 
expression 

_ (l/2)(3/2) ...(h + m + l/2)(l/2)(3/2) (g - m - 2 + 1/2) 
ym{P,ti)-y m (p + h)h\(p-l)\ 

where ip m is the value of the lowest formfactor. We considered the case p > m. One 
can prove that the case p < m leads to the same expression. Introducing the particle 
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momenta p — > —p, p < 0, and the hole momenta q = h the last formula takes the 
following form: 



1>m(p, ?) = i>: 



i r(-p-c) r(g + i + c ) 



(21) 



m 



(p-q)T(-p)T(l-c)T(q+l)T(cy 



where the parameter c equals 



1 

c = - + m. 



In the same way one can easily prove that for n particles and n holes the formfactor 
takes the form ip m F c (Pi,qi) where the function F c (pi,qi) is given by the equation (11) 
with the same parameter c = 1/2 + m. Again we observe that the expression for the total 
formfactor coincides with the prediction of the Luttinger liquid theory according to the 
operator equation (3). In fact the equation (3) gives c = y/£/2 + m/^ which coincides 
with our value of c at £ = 1. Let us stress once more that we have calculated the sets 
of the formactors corresponding to the m-th term in the correlator (5). Thus for the 
XX- spin chain we obtained the particle-hole formfactors for the low-lying excitations in 
agreement with the universal formulas predicted by the Luttinger liquid theory. Note 
that once the particle-hole formfactors for different m are obtained one can calculate the 
correlator using the formulas (9), (12) without mentioning the bosonization at all. Thus 
it is clear that the results of Section 3 for the correlator (prediction of the Luttinger 
liquid theory) is nothing else as the contribution of the complete set of the intermediate 
states corresponding to the low excitation energy. That means that in general the con- 
tributions of order ~ 1/L to the formfactors at low Pi, q,i ~ 1 can give the contributions 
to the correlator which are not suppressed by a power of L. Equivalently, in order to 
obtain the whole asymptotic expansion of the correlator one should take into account 
the contributions of the high-energy excited states. 

5. Exact expansion of the correlator. 

In this section we present the method to calculate the exact asymptotic expansion of 
the correlator up to the terms of arbitrary order. The correlator in the XX- spin chain 
G(x) which can be represented as a Toeplitz determinant [12]. In fact using the Jordan- 
Wigner transformation relating spin operators to the Fermi operators (af, a^) a+ = 
exp(i7r J2i< x n i) a t ) tlie correlation function G(x) can be represented as the following 
average for over the free-fermion ground state: 



where N(x) = J2i=i n i- Introducing the operators, which anticommute at different sites, 



G(x) = (0\a+e™ N Wa \0), 



A = 



Bi = aj - a,i 



AiBi = e 
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where rii = afdi - is the fermion occupation number, with the following correlators with 
respect to the free-fermion vacuum, 

(0|iVyo> = 2G (i-j) - Sij, (0|A^|0) = 0, (01^10) = 0, 

where Go(i — j) is the free-fermion Green function, one obtains the following expression 
for the bosonic correlator: 

G(x) = ^(0\B (A 1 B 1 )(A 2 B 2 ) . . . (A^B^A^O). 

Clearly in the thermodynamic limit the correlator is given by the determinant of the 
Toeplitz matrix M(i — j): 

G(x) = idety {M{i-j)) , i,j = l...x, 

where the matrix M{i — j) corresponds to the following generating function: f(x) = 
e ia; sign(7r/2 — \x\), —ir < x < n. Due to the form of the matrix M(i — j) (M(l) = for 
I- odd) the correlator G(x) can be represented in the following form [12]: 

G{x) = - (R N f , (x = 2N), G{x) = ^R N R N+1 , (x = 2N + l), (22) 

where R^- is the determinant of the N x N matrix of the following form: 

R N = det tJ ((-l) i -J'G (2i - 2j - 1)) , i,j = l,... N. 

Since Rjy is the Cauchy determinant in thermodynamic limit we obtain the following 
expression: 

The sum corresponding to the logarithm of the determinant R N was calculated in Ref.[15] 
in the context of the Ising model in the following form: 

1 oo (r\2k 1 \ JD 

H R N ) = ^A-- lnN + ^_- 2 J N ^ iy (24) 

where A = 2 1 / 12 e 3< >' ( ~ 1 * ) = 0.6450024... and B 2 k are the Bernoulli numbers. Using the 
expression (24) one can obtain from the equation (22) an arbitrary number of terms in 
the asymptotic expansion of the correlator (for example, see [16]). The first few terms 
are: 

«M - %^ ~ <-i)4i + 4i + (-i)4i - hh ~ 363 1 



y/jfx 1 / 2 ^ v ' Sx 2 128a; 4 v ' 8a; 4 64a; 6 v ' 1024a; 6 

1707 1 . lN:c 1985 1 
+ 32768^ + 1024^ + --^ (25) 
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One can see that only the first two terms in the equation (18) coincide with the exact 
result (25) i.e. only the first two terms in the expansion of the correlator correspond to 
the sum over the low-lying intermediate excited states. 

Conclusion. 

Using the expressions for the formfactors of local operators for the XX - quantum spin 
chain as a Cauchy determinants we obtained the contribution of the low-energy excited 
intermediate states to the asymptotic expansion of the correlator in the XX spin chain. 
We obtained the particle-hole formfactors corresponding to the higher order terms in the 
asymptotics of the correlator in agreement with the predictions of the Luttinger liquid 
theory. We have shown that in order to obtain the correct asymptotic expansion of the 
correlator in the XX- spin chain one should take into account the intermediate states 
with the high (of order of the band width) excitation energy. In the other words the 
contributions to the original formfactors which are suppressed by the power of 1/L at 
small quantum numbers of the particles and the holes, Pi,qi ~ 1, leads to the power-law 
contributions to the asymptotics of the correlator, which are not suppressed by a power 
of L. Equivalently, if various irrelevant operators are included into the equation (3), 
the sum over the corresponding low-energy formfactors (which are suppressed at low 
energies) can be saturated by the quantum numbers Pi, qi of order of L. Therefore the 
corresponding contributions cannot be described in the framework of the Luttinger liquid 
theory. We have seen that only the first two terms are in agreement with predictions 
of the Luttinger liquid theory. We expect the same results for the general case of the 
XXZ- spin chain. 
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Appendix. 

Let us calculate the matrix elements (8), (11) which correspond to the particle-hole 
formfactors. Let us begin with the matrix element (8) corresponding to the first branch 
of the Luttinger model [11]. We start with the matrix element for the operators in the 
coordinate space {0\a + (x)a(y)e Ba |0), B a = a(2Ti / L) J2 p >o Pi(p)/p- Using the formula 
Ae B = e s (J]n(l/ n [A-S] n )> one can commute the exponent e Ba to the left. Thus we 
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obtain the following equation: 

/ 1 _ P & \ a i / 1 - e ix ' \ a pw' 

(0\a+(x)a(y)e B «\0) = ^-^j (0\a+(x)a(y)\0) = - ^-^j 
where x' = 2nx/L, y' = 2ny/L. Then for the derivative we obtain 

(id x + i d y )(0\a + (x)a(y)e B «\0) = -a^ (l - e^J' 1 e« (l - e^')'^ . 

Calculating the Fourier transform for both sides of this equation we obtain for the matrix 
element (0|a+a p e Ba |0) the expression (8) for n — 1 with the factors 

/ + w=«r^"*" i)B ( i - ei9 )" ( " +i) - ?>o - 

27r dx „■„„. / Sn \a-\ 



where p and q are integers. Calculating the integrals and using the Wick's theorem we 
obtain the equation (8) with the functions f + (p), f~(o) presented in the text. 

The derivation of the equation (11) is similar. We consider the matrix element 
(0\c + (x)c(y)e Bc \0) with B c = c(2ir/L) J2 P <o P2{p)/p- Commuting e Bc to the left we 
obtain: 

/~\—p~w'\ c i f -\ — p-w' \° p-iy' 

(0\c+(x)c(y)e»°\0) = (j— (0\c + (x)c(y)\0) = ^ [j—,) ^ _ ^ ■ 

Calculation of the derivatives gives 

(id x + td y )(0\c + (x)c(y)e B c\0) = -c^ (l - e^J^e^ (l - e^)"^ . 

Calculating the Fourier transform we obtain for the matrix element (0|c+c p e Bc |0) the 
expression (11) for n — 1 with the factors 

/+(p> = r §i e " i " +i " j i 1 - e "")" t • p < "• 

r( q ) = c[' d ^(l-e-)-^\ 9 >0. 

Calculating the integrals we obtain the equation (11) with the functions f + (p), f~(q) 
presented in the text. 
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